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Abstract 

We summarize the current theoretical and experimental status of the spectral 
changes of vector mesons (p, uj, <ft) in nuclear medium. Various approaches includ- 
ing QCD sum rules, effective theory of hadrons and bag models show decreasing of the 
vector meson masses in nuclear matter. Possibility to detect the mass shift through 
lepton pairs in 7 — A, p — A and A — A reactions are also discussed. 
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1 Introduction 



At high temperature (T) and density (p), hadronic matter is expected to undergo 
a phase transition to the quark-gluon plasma PJ. The order parameter charac- 
terizing the transition is the chiral quark condensate (qq), the absolute value of 
which decreases as (T,p) increases. Numerical simulations of quantum chromo- 
dynamics (QCD) on the lattice are actively pursued to determine the precise 
nature of the transition at finite T 0] and various model calculations have been 
done to look for the observable signature of the phase transition ||. 

In this article, we will concentrate on one of the interesting critical phenomena 
associated with the QCD phase transition, namely the spectral change of hadrons, 
in particular the mass shift of light vector-mesons (p, to and 0) in nuclear matter 
at zero T. The vector mesons are unique in the sense that they decay into 
lepton pairs (e + e~ and p + p~) which can be detected experimentally without 
much disturbance by complicated hadronic interactions. 

In section 2, we will review the current knowledge of the quark condensate in 
medium. In section 3, various approaches to calculate the vector meson masses 
in nuclear matter are summarized. Section 4 and 5 are devoted to the detailed 
explanation of the mass shift of p, u and <fi in quantum hadrodynamics, an 
effective theory of mesons and baryons. Experimental possibilities to detect the 
spectral change are discussed in section 6. Concluding remarks are given in 
section 7. 



2 Quark condensates in nuclear matter 

The dynamical breaking of chiral symmetry has close similarity with the gauge- 
symmetry breaking in superconductors as has been first recognized by Nambu 



and Jona-Lasinio H. Table 1 shows a brief comparison between QCD and the 
BCS theory for "low temperature" superconductivity. 



QCD 


BCS 


q — q paring 


e T — e l paring 


chiral symmetry breaking 
(qq) ^ 


gauge symmetry breaking 
M ^ 


quark spectrum: E = ^Jp l + M' 2 
M (constituent quark- mass) ~ 350 MeV 


electron spectrum: E = ^Je(p) 2 + A 2 
A (BCS gap) ~ 0.01 eV 


Nambu-Goldstone (NG) boson = pion 


NG boson is absorbed by photon 



Table 1: Comparison between the dynamical breaking of chiral symmetry in 
QCD and the dynamical breaking of gauge symmetry in the BCS theory. 

The medium modification of the quark condensate has been calculated by 
various methods (lattice QCD, chiral perturbation theory, Nambu- Jona-Lasinio 
model etc). See a review and also ||. By these studies, it turned out that 
there is one noticeable difference between the behavior of (qq) at finite T (with 
p = 0) and that at finite p (with T = 0): In the former case, the significant 
change of the condensate can be seen only near the critical point T ~ T c J7]]. 
On the other hand, in the latter case, O(30%) change of (qq) could be seen even 
in normal nuclear-matter density. This observation is based on the following 
formula in the fermi-gas approximation (independent particle approximation) || 
{uu) 4S^jv fP F d 3 p M N 

(uu) ~ f 2 mlJ (2tt)3£(p) ' 1 j 
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Here m^m^) is the nucleoli (pion) mass, f n is the pion decay constant, T, n 



N 



(45 ± 10)MeV is the ttN sigma term, and E(p) = \jp 2 + M%. (■) and (-) denote 
the expectation value in nuclear matter and that in the vacuum respectively. 
The integration for the nucleon momentum p should be taken from to the fermi 
momentum p F . At normal nuclear matter density (p = po = 0.17/fm 3 ), the above 
formula gives (34±8)% reduction of the chiral condensate from the vacuum value. 
In Fig.l, (uu) I (uu)o as well as the strangeness condensate (ss) / (ss) are shown in 
the linear density approximation ||, where the uncertainty of E^-jy is considered. 
Estimates taking into account the fermi motion and the nuclear correlatons show 
that these corrections at p = p are less than the above uncertainty . 



Unfortunately, the condensate itself is not a direct observable and one has to 
look for physical quantities which are measurable and simultaneously sensitive to 
the change of the condensate. The masses of light vector-mesons are the leading 
candidates of such quantities. 

Fig.l 



3 Vector mesons in nuclear matter — overview 



Let's consider p, ui and mesons propagating inside the nuclear matter. Adopting 
the same fermi-gas approximation with and taking the vector meson at rest 
(q = 0), one can generally write the mass-squared shift as 

r 2 *2 2 „ f PF d3 P M N n ( X / Q X 

6m v =m v -m v = Aj __/^(p), (2) 

where /vjv(p) denotes the vector- meson (V) - nucleon (N) forward scattering 
amplitude in the relativistic normalization, and my (my) denotes the vector me- 



son mass in nuclear matter (vacuum). Here, we took spin-isospin average for 
the nucleon states in fvN- If one can calculate /yjv(p) reasonably well in the 
range < p < p F = 270 MeV (or 1709 MeV < y/s < 1726 MeV in terms of the 
V — N invariant mass), one can predict the mass shift. Unfortunately, this is a 
difficult task: /yjv(p) is n °t a constant in the above range since there are at least 
two s-channel resonances iV(1710), iV(1720) in the above interval and two nearby 



resonances iV(1700) and A(1700). They all couple to the p — N system |l ij and 
give variation of /vjv(p) as a function of p in principle. From this reason, one 
should develop other methods to estimate 5my without refering to the detailed 
form of /yjv(p)- We will briefly review two of such approaches in the following 
subsections, namely the QCD sum rules and effective theories of hadron. 

3.1 QCD sum rules 

This subsection is partly based on the work in ref . [|T2"| , |13[ . The QCD sum rules 
(QSR) for vector mesons in nuclear matter were first developed by Hatsuda and 
Lee []12[. In their approach, one starts with the retarded current correlation 



function in nuclear matter, 

IV(u;,q) = i J dW^RJ^aOJ^O)) , (3) 

where g M = (u, q) and RJ M (x) J u (0) = 9(x°)[J^(x), J v (0)] with the source currents 
J M defined as J? ,UJ = ^(uj^u =F dj^d) (— (+) is for the p (cj)-meson) and J* = 
s7 M s. Although there are two independent invariants in medium (transverse 
and longitudinal polarization), they coincide in the limit q — >• and reduce to 
n^/f- 3a; 2 ) = II. II satisfies the following dispersion relation, 

1 r°° ImYl(u) 
Ren(c; 2 ) = -P / du 2 , K i + (subtraction). (4) 

71 JO U — UJ 2 
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In QSR, the spectral density Imll is modeled with several phenomenological pa- 
rameters, while Rell is calculated using the operator product expansion (OPE). 
The phenomenological parameters are then extracted by matching the left and 
right hand side of (f§) in the asymptotic region u 2 — > — oo. The density depen- 
dence in the OPE side is solely determined by the density dependent condensates 
which are evaluated from low energy theorems or from the parton distribution of 



the nucleon [fl2 | . 



In the medium, we have three kinds of structure in the spectral density: the 
resonance poles, the continuum and the Landau damping contribution. For q — > 
0, the last contribution is calculable exactly and behaves like a pole at uo 2 = 



12, 14]. In total, the hadronic spectral function looks as 



87rImn(u>(T) = 5{u 2 )p sc + F*6(u 2 - m 2 v *) + (1 + — )9(u 2 - S* ) (5) 



7T 



= Phad.{u 2 ), 



with p sc = 2n 2 p/^Jp 2 F + Mjf ~ 2n 2 p/M N . m v , F* and Sq are the three phe- 
nomenological parameters in nuclear matter to be determined by the sum rules. 

Matching the OPE side and the phenomenological side via the dispersion rela- 
tion in the asymptotic region uj 2 — > — oo, we can relate the resonance parameters 
to the density dependent condensates. There are two major procedures for this 
matching, namely the Borel sum rules (BSR) |15| and the finite energy sum rules 
(FESR) jl6| , which can be summarized as 

/•oo 

/ dsW(s) \phadis) -p OPE {s)} = 0, (6) 
Jo 

s n 9(S -s) (FESR), 

e s/M 2 

Here the spectral function phad.i. 3 ) stands for eq.(|5|). p OPE (s) is a hypothetical 
imaginary part of II obtained from OPE. 
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W(s) 



(BSR). 



To make quantitative analyses of spectral parameters, the stability analysis 
based on the Borel transform is more suitable than FESR. Since the Borel mass 
M is a fictitious parameter introduced in the sum rule, the physical quantities 
should be insensitive to the change of M within a Borel interval M min < M < 
M max ; namely the principle of minimum sensitivity (PMS) [I?]] is used. One 
can accomplish this insensitivity by choosing Sq suitably at given density. In 
Fig. 2, the Borel curves for the p(u) meson for three different values of baryon 
density are shown with Sq chosen to make the Borel curve as flat as possible in 
the interval 0.41GeV 2 < M 2 < 1.30GeV 2 . The upper (lower) bound of the Borel 
interval is determined so that the power (continuum) correction after the Borel 
tranform does not exceed 30 % of the lowest order term in OPE. 

Fig.2 

By making a linear fit of the result, one obtains [O, 13 



m 



p,LU 



P 



and 



1 - (0.16 ±0.06) — , (7) 

m p ,u> Po 

\hr = 1- (0.15 ±0.05) A (8) 
V Po 

— = 1 - (0.24 ±0.07) — , (9) 
F Po 



^ = 1 - (0.15 ±0.05) y — , (10) 



where y is the OZI breaking parameter in QCD defined as y — 2(ss) N / (uu + dd) N 
with (-)jv being the nucleon matrix element, y takes the value 0.1 — 0.2 fl~2[ 



The decrease in eqs. ([7|,|T0|) is dictated by the density dependent condensates 
(qq), ({qq) 2 ) and (q , ~f fl D u q). The errors in the above formulas are originating 
from the uncertainties of the density dependence of the these condensates. The 



contribution of the quark-gluon mixed operator with twist 4, |H1 which may 



possibly weaken the mass shift, is neglected in the above. Shown in Fig. 3 is the 
mass shift given in eqs. ([7|,[10]) with possible theoretical uncertainties. 

Fig.3 

Some sophistications of the QSR analyses by Hatsuda and Lee have been done 
later by several authors. 

(i) Asakawa and Ko have introduced a more realistic spectral function than @ 
by taking into account the width of the p-meson and the effect ofn — N — A — p 
dynamics fl8f . By doing the similar QSR analysis as above, they found that the 
negative mass shift occurs also in this realistic case. 

(ii) Monte Cairo based error analysis was applied to the Borel sum rule by Jin 



and Leinweber [|19[ instead of the Borel stability or PMS analysis employed in 
P, g. They found m*Jm Pfi) = 1 - (0.22 ± 0.08) (p/p ) and m;/m = 1 - 
(0.01 ± 0.01) (p/po), which are consistent with eqs. (fFHTUD within the error bars, 
(hi) Koike analysed an effective scattering amplitude Jvn defined as 8m 2 v = Jvn'P 
using the QSR in the vacuum |2TJ. Although his original calculation predicting 



fv n > is in error as was pointed out in ref. [13, [19[, revised calculation gives a 



consistent result with eqs. (|7| JT0|) within the error bars |2lJ. Note here that /i 



IVN 



does not have direct relation to the scattering length at zero momentum /Vtv(0). 
3.2 Effective theories 

There have been many attempts so far to calculate the spectral change of the 
vector mesons using effective theories of QCD. The first attempt by Chin [2l| 
using the quantum hadrodynamics (QHD) shows increasing tu-meson mass in 
medium due to a process analogous to the Compton scattering; 

u + N -> u + N. (11) 
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For the p-meson, similar but more sophisticated calculations taking into account 
A-resonance and in-medium pion show a slight increase of the p-meson mass 
[E3|. In these calculations, only the polarization of the Fermi sea (the particle- 
hole excitations) was considered. Also their predictions are different from the 
general assertion by Brown and Rho claiming that all the hadron masses except 



for pion should decrease [24 



On the other hand, Saito, Maruyama and Soutome [53] and Kurasawa and 



Suzuki p6| have realized that the mass of the cu-meson is affected substantially 
by the vacuum polarization of the nucleon in medium 

io -> N*N* -> u, (12) 

where N* is the nucleon in nuclear matter which has smaller effective mass than 
that in the vacuum. They show that the vacuum polarization dominates over the 
Fermi-sea polarization in QHD and leads decreasing vector meson mass. This 



conclusion was later confirmed by several groups [27, ^8], |29| and was generalized 



for the p and (j) mesons by the present authors |31| which will be discussed in 
more details in section 4 and 5. Jaminon and Ripka has also reached a similar 
conclusion by using a model of vector mesons coupled to constituent quarks |32|] . 



Saito and Thomas have examined a rather different but comprehensive model 
(bag model combined with QHD) and found decreasing vector-meson masses 
iHi m * P ,u/ m p,u ~ 1 — 0.09(p/po)- The spectral shift of the quarks inside the bag 



induced by the existence of nuclear medium plays a key role in this approach. 

Basic idea common in the approaches predicting the decreasing mass may 
be summarized as follows. In nuclear matter, scalar (a) and vector (u) mean- 
fields are induced by the nucleon sources. These mean-fields give back-reactions 
to the nucleon propagation in nuclear matter and modify its self-energy. This 
is an origin of the effective nucleon mass < in the relativistic models 
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for nuclear matter. The same mean-fields should also affect the propagation of 
vector mesons in nuclear medium. In QSR, the quark condensates act on the 
quark propagator as density dependent mean-fields. In QHD, the coupling of 
the mean-field with the vector mesons are taken into account through the short 
distant nucleon loop with the effective mass M^. In the bag- model, the mean 
fields outside the bag acts on quarks confined in the bag and change their energy 
spectrum. 

Let us show here that one can understand the negative mass shift of the vector 
mesons in a simple and intuitive way in the context of QHD. More quantitative 
discussion will be given in the later section. After renormalizing infinities in the 
vacuum loop, the density-dependent part of the Dirac-sea polarization to the 
vector-meson propagator is approximately written as 

D{q) ~ Z- V _ ml = qi-Zmy (13) 

where Z being the finite wave-function renormalization constant in medium. The 
pole position is thus obtained as m v = y/Zmy- Because of the current conser- 
vation, only the wave function part of the propagator is modified in medium. 
Since the effective mass of the nucleon decreases in medium (M^/Mjv < 1), 
physical vector mesons have more probability to be in virtual baryon— anti- 
baryon pairs compared to that in the vacuum. This means Z < 1, which leads 



to m* v /m v = Z < 1 [2JL B0| 



M* N /M N < 1 -> Z < 1 -> m* v /m v = Z<1 . (14) 



4 p and uj mesons in quantum hadrodynamics 
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4.1 Nucleon at finite density 

Before discussing the vector meson masses in medium, let's make a brief review 
of the effective nucleon mass at finite density in the quantum hadrodynamics 
(QHD) on the basis of ref. [34, |22|, |35|. The lagrangian of QHD is written as 



C = $[-Yn(id» - g v V») - {M N - g 8 S))i/> 

+ ^(d,Sd»S -m 2 s S 2 ) -\f, u F^ + ~mlV t y» + C C r , (15) 



with 



F lw = d i y v -d v V lt , (16) 

where ip(x), S(x) and V(x) are nucleon, a meson and u meson field, respectively. 
u{uj) field is coupled to nucleon current with the coupling constant g s (g v )- C-ct 
is the counter term added to the original lagrangian in order to get the finite 
physical quantities. 

Let's define the nucleon propagator in uniform nuclear matter, 

iG a p(x,y) = (Tip a (x)^(y)) . (17) 

The free nucleon propagator in momentum space is expressed as 

GV ) = (7 ■ P + M N ) [ - 2 _ ^ - - + JL. S(P0 - EMM, - |p|)} , 

= G F (p) + G° D (p), (18) 



where E{p) = yp 2 + Mff and pp is the fermi momentum. In the relativistic 
Hartree approximation (RHA) |35|, |22j , the full nucleon propagator reads 



G*(p) = G°(p) + G°(p)J: H G H (p), (19) 

where H denotes the Hartree approximation. A schematic diagram for G H (p) is 
given in Fig. 4. 
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= K- 


7^ 




= 4 


/" A 


mf 


y (2tt) 4 




= 4 








y (2tt) 4 



Fig.4 

E ff in Fig.4 can be written as 

(20) 

Tr[G H (p)} } (21) 
Ti[YG H (p)}, (22) 

where S S (S^) is scaler (vector) self-energy. Formal solution of eq.(plj) reads 

[G H {p)\- 1 = 1 -p-M N -Z H 

= 7 -(p + S i; )-(M A r + S s ), (23) 

or equivalently 

G"{p) = ( 7 ■ p + m [ ¥ _^ + ie + j^Sfr - E-<p))0<p F - |p|)} , 
= G H F {p) + G H D {p), (24) 



where = ^/p 2 + M^ 2 ,p = p + £„ and = M N + S s . Eq. (gD implies 

that the interacting propagator can be separated into two parts: Gp(p) which 
has the same form as the free nucleon propagator with the effective nucleon mass 
Mjy, and G%{p) which depends explicitly on the fermi momentum in nuclear 
matter. The advantages of this separation will be discussed later. 

Next we turn to discuss self-energy S H in which E s is related to the nucleon 
effective mass. For the vector self-energy, insertion of eq.(p4"|) into eq.(p2"D gives, 



m 2 v J (2tt 4 ) p 2 - M$ + ie m 2 

where Pb = ^P 3 f * s ca hed baryon density with 7 = 4 for nuclear matter. By 
shifting the integration variables from p to p, the first term on the r.h.s. vanishes, 
namely has no divergent term. In the case of scalar self-energy, however, one 
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.gl f d*p 8M* N _ gl 7 fp r M* N 



gets 



'mji (27r) 4 p 2 -M^ 2 + ie m 2 (2vr) 3 7 K (p 2 + M^ 2 ) 1 ^ 

where the first term on the r.h.s. is divergent. To single out this divergence, we 
extract the infinite part using dimensional regularization: 

^nf = _^£(2_|) 

= \ 2 2 {M% + 3M 2 / E S + 3MArSg + X^), (27) 

where £* denotes the infinite parts and the diagrammatic illustrations are de- 
picted in Fig. 5. 

Fig. 5 

These divergent terms are removed by defining the following Cct] 

C CT = ai S + ^ + p3 + ^S\ (28) 
which yields new contributions, E^ T , to the scalar self-energy: 

Sf r = i;A( Z T)(^ i ) B «^i- (29) 

The corresponding Feynman diagrams are depicted in Fig.6. 

Fig. 6 

a\ — «4 are determined so as to cancel precisely the divergent parts in (p7|): 



«n = hx-^H" - 1)! / 4^n G F(pn (so) 



Note that the finite parts of ai_4 are also fixed in the above conditions 
], which may or may not be justified and must be checked using experimen- 
tal/empirical inputs. Examining this point is not a main theme of this article, 
and we simply take the above procedure to calculate the density dependent part 
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of 



4^ + 4h» M^Mv) 



-M 2 N (M* N - M N ) - -M N (M* N - M N ) 



M* N - M N , 



Y {M * N 



M 



Nj 



(31) 



where 



M 



Ps 



N 



71- 



p F E*(p F )-M* N 2 log 



p F + E*{p F ) 



M* N 



(32) 



This expresses the effective nucleon mass including the vacuum polarization. 
The density dependence of the effective nucleon mass is shown in Fig.9. g s ,g v 
and m s are chosen to satisfy the saturation density for nuclear matter (—15.75 
MeV) at nuclear matter density p F = 1.30 fm~ : g s = 7.37, g v = 10.1 and 



m s = 458MeV p9|. Hence we find that effective nucleon mass has the reduction, 



M^/Mtv = 0.730, at nuclear matter density. As will be shown in section. |4T^ , the 
reduction plays a crucial role in studying vector meson masses at finite density. 



4.2 p and u mesons at finite density 

This subsection is partly based on the work in ref. |3(| 



4.2.1 Effective lagrangian 



Let's start with an interaction lagrangian of p, uj with the nucleon: 



Lint gv 



Vf, v = {p,u} 



(33) 



where a runs from through 3, Vq (Vis) corresponds to the u (p) field, r a is the 
isospin matrix with t°=1, and is the nucleon mass. The numerical values of 
the coupling constants (g v , k v ) will be given in sec. 4.2.4. 
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In nuclear medium, meson propagator is defined as 

jr( I)y ) = {rv(i)r(y)). (34) 

By using the free vector-meson propagator 

D% u (q) = 9flu + 1 - - , (35) 

" ' q 2 — m? + ie A (q 2 — m?/\ + ie)(q 2 — m? + it) ' 



with A being a gauge parameter in the Steukelberg formalism |41| , the full prop- 
agator can be written as 

D""{q) = D^(q) + D^(q)U Xa (q) D~ (q) , (36) 

with the self-energy H^ u . See Fig. 7. a. 

Fig.7 

In the one-loop level, the density dependent part of the self-energy comes only 
from the nucleon-loop (Fig.7.b): 

K (?) = -7^)4 / d'kTv[TlG H {k + q)t b v G H {k)) , (37) 
where (a, b) are the isospin indices and we have used G H (p) defined in section 



4.1. For the vertices, we make use of 



^ = 9vh,T a -^-*» x iq X T% ft = 9v [l„T h + Jj^igV] • (38) 

where the relative sign of the tensor part in and that in T v is opposite to that 
in ref. [^6|. The self interaction of the p meson gives density dependence only 



from two or higher loops. The coupling of p with in-medium pions analyzed in 



23|| is also the higher loop effect and will not be considered in this paper. 



Ii^ u in fl37|) is composed of two parts 11^ = H ^ + n£,: the first term corre- 



sponds to the fluctuation of the Dirac sea of the nucleons with mass Mjy-, while the 
second term gives the fluctuation of the Fermi sea + the Pauli blocking. The ad- 
vantage of taking this separation is that 11^ and ITd satisfy current conservation 

15 



separately, i.e., pTI°f (p) = pTI^, = PBfl . If one adopts a separation without 



current conservation, the corresponding self-energy induces spurious results |38|] . 
11^ generally has divergences to be subtracted. We will show our subtraction 
procedure in section |4.2.2| and define as the subtracted polarization. 



The vector meson propagator in the medium has a general form 

Dflu = hrrr + ~ frrr > ( 39 ) 

qz _ m z _|_ gz _ m z _|_ ^ 

where we have suppressed isospin indices (a, b) and m denotes the p or u mass in 
the vacuum. P^ v {P^) is the projection operator to the transverse (longitudinal) 
direction to q: 

= g^(9v + m/ci 2 )g u >, P^ = ^e v (40) 

with e " = -^(|qUoq/|q|)- (41) 

Vq 2 

Ht,l is related to U^ v as 

Hz, = -(g 2 /q 2 )n 00 , n T = (n{ + (g 2 / q 2 )n o)/2. (42) 

In this paper we will focus on Ht,l in the time like region with q = where IT^ 
is exactly equal to Hp. To obtain (p9[), we adopt (|35| ) with A — > oo as a free 
propagator of the massive vector mesons. 

4.2.2 Dirac sea effect and subtraction procedure 



The interaction ( p3|) is not a renormalizable one in the conventional sense. This 
implies that the model contains infinite series of the higher dimensional operators 
which play a role to cancel the divergences emerging from the loops of the lower 



dimensional operators [42, 43]. Instead of developing a systematic subtraction 



procedure, we will take a phenomeno logical way to extract U^ u from H°J^. First of 
all, we will normalize the propagator at zero density as l/{q 2 — m 2 ), i.e. subtract 
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away both the divergent and finite parts from IT^. This corresponds to a set of 
the renormalization conditions at zero density, 
d n U F (q 2 ) 



(n = 0,1,2,- • -,oo). (43) 



d(q 2 ) n 

A straightforward generalization of the above conditions to finite density reads 
d n U F (q 2 



= (n = 0,1,2,- • -,oo). (44) 

M* ->M N ,q 2 =m 2 



d(q 2 ) n 

This together with a requirement that the higher dimensional counter terms are 
the local polynomials allow one to single out the density dependent part of 
uniquely. 

Our renormalization condition is different from the simple scheme U F U = 

jOF ttOF 



— 11^ \m* n =m n - This simple scheme cannot remove the divergences such as 
M^r/e and (M^) 2 /e (e — > 0) induced by the pNN tensor coupling. These diver- 
gences require the counter terms such as S n F flu F^ lu (S being the scalar field and 

being the field-strength for the vector fields) which are allowed in our scheme 
but not in the simple scheme. Although our procedure is physically plausible, 
it is still "a" way to subtract the divergences among many other possibilities. 



For the to meson, our procedure is equivalent to that in PS). For m*, we have 



checked that the different subtraction procedures in [26, [29| does not cause more 
than 3% differences from ours at p = Pq. This is shown in Fig. 8. 

Fig.8 

Using the dimensional reguralization and the above subtraction procedure, 
one obtains the following TL F for the p and to mesons. 

K F = ^(g^/g 2 -g„)(u F + n^ + nf) 

ttF 9v 2 f 1 , m \ M* N 2 - q 2 x{l - x)\ 

U F = ^q 2 dxx(l-x)\og\-% r . (45) 



7r 2 Jo { Mfj- — q 2 x(l — x) 

( g 2 v K v ^M* N q 2 /-i Jml _\M* N 2 -g 2 x{l-x) 



<> - iw'-fji^ U-w-.) !' (46) 
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^ = ( 2M } 2W0 ^ {M " +gX(1 ~ X)}1 ° g l M^- g Ml-x) | G47) 
where the integrals take analytic forms: 



A* f2M*J \ , t „ A (2Ml \ i , . . 

+ T f + 1 j tan-^l/A*) - - + lj tan-^l/A), (48) 

Jo I Mx - q 2 x(l - x) J 

= 2\og(M* N /M N ) + 2A*ttm- 1 (l/A*) - 2Ata,n~ 1 (l/A), (49) 

where A* = - 1|V * , A= - 1|Y and 4Mff > q 2 is assumed. 



4.2.3 Fermi sea effect 

Another contribution to IL^ comes from the Fermi sea and the Pauli blocking 
effects, which is expressed as 

K°^) = -7^)4 / d'mVlG F {k + q)T b u G D (k) + (F «-> D)] (50) 
= ^K + nJ + nf^. (51) 



Substituting eq.(|24]) into eq. (|50|) , one obtains the following form: 

(52) 

4M^g 2 J (g), (53) 















(nf 




/ \ 2 < 



-n^(g) + Q / ,J(4M^ 2 + g 2 )/ (g)+ ,K 



M 

where = {q^q u /q 2 — g^u) and p s has been defined in fl3"2]), and 
VWJ ~ W E*{k) g 4 -4(fc-g) 2 ' 



,(54) 
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x {{k-qf 



q 2 



+ q z 



k-q 



i n (q) = J n (q) + Jn(-q), 

j ,s _ -2 f d 3 k e(p F -k) (2E*(k) + q ) n 



I, k ' q 



k =E*(k) ■ 



, (55) 
(56) 
(57) 



(2tt) 3 J E*(k) q 2 + 2q-k + ie 
Here we have shown only the real parts of I1 D relevant to study the pole position 
at q = 0. 

In the limit of q = 0, Io(q) and ITr^g) defined from Yi^q) read 



n T (g ,q = 0) = n L (g ,q = 0) 

E*(p F )+p F 



/o(go,q = 0) = ^|21og 



M 



N 



M$ + % | A* tan" 1 



Pf 



E*(p F )A* 



p F + E*{p F ) 



M 



N 



- 2A*tan~ 1 



Pf 



E*(p F )A* 



(58) 



(59) 



where A* 



4 < _ 1 

4 



1/2 



and 4M* 2 > q 2 , is assumed. 



The mixing of the to meson and the a meson does not contribute to the 
propagator as far as q = and q^ > 0. This is because the meson self-energy 
relevant for the u — a mixing vanishes exactly in the time-like region with q = 



4.2.4 Coupling constant 

Since the uNN tensor coupling is generally small (e.g. k w = 0.12 in the vector 
dominance model), we take (g u , k u ) = (10.1, 0.0) as a typical strength. For the p 
meson, we adopt the following two sets. 





set I 


set II 


9 P 


2.63 


2.72 




6.0 


3.7 
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Table 2: Two sets of pNN coupling constants. g p (k p ) denotes the vector (tensor) 
coupling. 



Set I is obtained from the N — N forward dispersion relation [g4| . The Bonn 
potential of the N — N force gives similar values with this set. Set II is obtained 
by the vector-meson dominance together with the p universality E|. A major 



difference between the two sets is the strength of the pNN tensor coupling. (See 



ref.pO| for the detailed discussion on the vector- meson coupling constants.) In 



our calculations, vertex form factors are not taken into account for simplicity. 
4.2.5 Numerical Results 

As mentioned before, we will focus on the time like region and consider the inverse 
propagator 

Dj\q , q = 0) = q 2 - m 2 + Tl°(q , q = 0) + Il£(gg) . (60) 

Let us define two kinds of masses m* e (real mass) and m* nv (invariant mass). 
m* e is defined as a lowest zero of D^ 1 (q ,0). It is the quantity to be compared 
with that in the QCD sum rules. The invariant mass m* nv is defined as a lowest 
zero of D^ 1 with neglected, in which case D^ 1 is a function of q 2 only. m* nv 
here contains only the fluctuation of the Dirac sea by definition. 

In Fig. 9, the effective masses of to are shown together with M^/M^. The 
dashed line denotes m* nv /m. One sees that Both m* e and m* nv decrease at finite 
density, e.g. m* e /m ~ 0.8 at p = p . The behavior of m* einv confirms the 



importance of the Dirac sea polarization found in previous studies J37 , |26 



28], |29], |3(J. m* e < m re is caused by the fact that the physical uo is more dressed 



by the AW pairs in the medium since < Mn as we have discussed in sec. 3. 2 
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Fig.9 

In Fig. 10 and Fig. 11, we have shown the effective masses of the p meson with 
the parameter set I and set II, respectively. The strong pNN tensor coupling 
plays a dominant role and gives m* e /m ~ 0.6 — 0.7 at p — Pq. The polarization of 
the Dirac sea is again the most important ingredient and the suppression of the 
wave-function renormalization factor Z is the main reason for the mass reduction. 

Fig.10 and Fig.ll 

It is in order here to make some remarks: The reduction of m* e /m is consis- 
tent with that in the QCD sum rules (eq.(7)) for the u meson, and even larger 
reduction is observed for the p-meson in QHD. One should, however, notice that 
we neglected vertex form factors for the pNN and ujNN couplings. Such form 
factors will generally attenuate the magnitude of the mass shift of p and u. From 
Fig. 9-11, one also observes considerable non-linearity of function of 

density, which is contrast to the linear dependence in eq.(7). 



5 meson in generalized QHD 



This section is partly based on the work in ref . |31| . The 0-meson is a ss resonance 
in J p = 1~ channel with a narrow width (m^ = 1019.4 MeV and = 4.4 MeV). 
It is a suitable probe of the partial restoration of chiral symmetry in hot/dense 



hadronic matter together with the p and to mesons [46, 12, |4T| . In this section, 
we start with an effective hadronic model where the 0-meson couples to nucleon 
and hyperons (B = N, A , S 1 * 1 , E°) with the vector coupling, 



L in t = Y,9*B B-ypB (61) 

B 

g^B is the 0-baryon coupling constant listed in Table 3. 
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Baryons 


9aB 


9uiB 


9<t>B 


AT 


8.7 


10. o 


4.2 


A 


5.2 


6.9 


4.9 


S 


5.2 


6.9 


4.9 



Table 3: g^, g W B and g^e denote cx-i? scalar coupling, u-B vector coupling, and 
(f)-B vector coupling, respectively. g a B and g^B are taken from BE . 



Some remarks are in order here: 

(i) cf) — A and (f> — £ couplings do not break the OZI rule, while the <p — N coupling 
is OZI violating. 

(ii) H is neglected, since its effect to the self-energy is doubly suppressed by 
the mass of S and by the OZI violation in — H coupling. 

(hi) If one relies on the quark counting rule |[49|| , the 0-hyperon couplings are 
related to the tu-hyperon couplings as ^a(^e) = fcA( u s)/\/2 with s^a^e) being 
determined by the fit of the hypernuclear levels This is assumed in Table 3. 
(iv) OZI violating 0-nucleon coupling is not known experimentally. A study of 
the electromagnetic form-factors of the nucleon, however, yields an upper bound 



of its strength |50fl : g^/g^N < 0.4. 

For the non-strange nuclear matter, effects of the hyperons to the 0-meson 
self-energy arise only through hyperon— anti-hyperon loops, while the nucleon 
contribution to the self energy has both N — N loop and the Compton-type 
process. The one-loop self energy from hyperon and nucleon contributions reads 



KM) 



(2tt)< 



d'k g\ B Tr[ 7 ,G H (k + q) lu G H (k)} 



(62) 
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where G H (k) denotes baryon propagator in nuclear matter given in (^3) with 
M* N replaced by M* B . 

Although one can calculate the density dependence of Mg within QHD as in 
sec. 3.1, we take here the following simple ansatz to illustrate the essential feature 
of the correlation between M* B and m^: 

9aA{aH)/9aN = (M A(E) - M^) / (M N - M* N ), (63) 
M* N /M N ~ l-0.15(p/p ). (64) 

Eq . (|63"D is an universal relation in QHD |48] and Eq.(^) is a simple parametriza- 



tion of the effective nucleon mass which sometimes used in the literatures for 
p < 2p (see, e.g. ^8|). In Fig. 12, effective masses of N, A and £ parametrized 
by eqs. ([)3|,[>4]) are shown as a function of baryon density. 

Fig.12 

The </>-meson mass in medium is obtained as a solution of the dispersion 
relation 



io — m. 



+ £n B (^,0) = 0, (65) 

B 

where n B (u;,0) = — Yl B i (u, 0)/3u> 2 , and is the 0-meson mass in the vac- 
uum. IT^f is assumed to be a renormalized self-energy with the same subtraction 
procedure in the previous section. 

In Fig. 13, the ratio m^/m^ calculated with only the hyperon-loops is shown by 
the solid line, decreases by 6% at p = p . Note that the Okubo-Zweig-Iizuka 
(OZI) rule is preserved for 0-hyperon vertices, while it is violated in the self- 
energy n#. This is because the self-energy represents interaction of (ss pair) 
with non-strange nuclear matter. Similar phenomena are known in two-step decay 
processes such as (f> — » KK — > pn, f — > KK — > nn, and J/ip — > DD — > pn, 
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where each vertex preserves the OZI rule while the whole amplitude does violate 



the rule [5T 



Fig.13, and Fig.14 

The solid line in Fig.14 shows the ratio m^/m^ calculated with only the 
nucleon-loops. Since g^N is not known experimentally, the result in this case has 
much uncertainty compared to the hyperon case. We have used g^N/g^N = 0.32 
in Fig.14 which is close to the upper bound given in Table 3: thus the resultant 
decrease of in Fig.14 should be considered as an upper limit. The negative 
mass shift in Fig.13 and Fig.14 is a direct consequence of the current conservation 
(d^Bj^B) = 0) and M* B fM B < 1 as discussed in sec. 3. 2. This mechanism is 
quite general and does not depend on the details of the interaction and on the 
virtual particles running in the loop. 

To see the effect of the ultraviolet cutoff on the finite part of the loop integral 
in (HD, let us define n£ M (a>, 0; A^) = II^(u;, 0; A cut ) - W^(uj, 0; A cut ) | p=0 and 
use this in fl55p. We take covariant cutoff for A cut for simplicity. When A cut — > oo, 
0; A cut ) reduces to the renormalied IIb(u;,0). The dashed lines in Fig.13 
and Fig.14 are the results of such calculations for three cases, A cut = 1, 2, 10 
GeV. Although the cutoff dependence is not negligible, the qualitative picture 
we draw in the above is not affected. 

We have considered only the nucleon and hyperon loops in the </> self-energy. 
Another possible contribution is the kaon-loop in medium. This was studied by 



Ko et al. [47| who found that the kaon-loop has a tendency to decrease m*. at low 
densities provided that the effective kaon mass m* K = (m* K ^ +m* K+ )/2 decreases 
in medium. However, whether fh* K really decreases in nuclear matter or not is 



still a controversial issue (see e.g. @). 
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6 Experiments 



How one can detect the spectral change of vector mesons in experiments? One 
of the promising ideas is to use heavy nuclei and produce vector mesons in 7 — A 
or p — A reactions. Suppose one could create a vector meson at the center of a 
heavy nucleous. (It does not matter whether it is created at the nuclear surface 
or at the center as far as the produced vector mesons run through the nucleous 
before the hadronic decay). It is easy to see that the number of lepton pairs 
decaying inside the nucleous N iri {l + l~) and that outside the nucleous N out (l + l~) 
are related as 

N out (l+l-) ~ e-^ t R ' [ > 

where T to t denotes the total width of vector mesons ((1.3frn) -1 , (23fm) _1 and 
(45ml)- 1 for p, u and <f), respectively) and R being the nuclear radius. Eq.flBSf) 
shows that even the <fi meson has considerable fraction of N in /N out if the target 
nucleous is big enough. 

There exist already some proposals to look for the mass shift of vector mesons 



in nuclear medium |p3| . One is by Shimizu et al. who propose an experiment 
to create p and u in heavy nuclei using coherent photon - nucleus reaction and 
subsequently detect the lepton pairs from p and uj. Enyo et al. propose to 
create meson in heavy nuclei using the proton-nucleus reaction and to measure 
kaon pairs as well as the lepton pairs. By doing this, one can study not only 
the mass shift but also the change of the leptonic vs hadronic branching ratio 
r = T(0 — > e + e~)/r(0 — > K + K~). Since is very close to 2m^- in the vacuum, 
any modification of the 0-mass or the .fT-mass changes the ratio r substantially as 
a function of mass number of the target nucleous. Similar kinds of experiments 
are also planned at CEBAF and GSI. 
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There are also on-going heavy ion experiments at SPS (CERN) and AGS 
(BNL) where high density matter is likely to be formed. In particular, CERES/NA45 
and HELIOS-3 at CERN reported enhancement of the lepton pairs below the p 



resonance 



f>H |55|1 , which may not be explained by the conventional sources of 



lepton pairs. E859 at BNL- AGS reported a possible spectral change of the 0-peak 



in K + K~ spectrum ||56|| . If these effects are real, low mass enhancement of the 
lepton pair spectrum expected by the mass shift of the vector mesons could be a 
possible explanation |)7] (see also ]5l| for another explanation). In nuclear col- 
lisions at higher energies (RHIC and LHC), hot hadronic matter or possibly the 
quark-gluon-plasma with low baryon density are expected to be formed. In such 



cases, double 0-peak structure proposed by Asakawa and Ko [47] as well as the 



spectral change of p and uo |59| will be a distinct signal of the chiral restoration 
in QCD. 



7 Concluding remarks 

The spectral change of the elementary excitations in medium is an exciting new 
possibility in QCD. By studying such phenomenon, one can learn the structure 
of the hadrons and the QCD ground state at finite (T, p) simultaneously. Theo- 
retical approaches such as the QCD sum rules and the hadronic effective theories 
predict that the light vector mesons (p, u> and 0) are sensitive to the partial 
restoration of chiral symmetry in hot/dense medium. These mesons are good 
probes experimentally, since they decay into lepton pairs which penetrate the 
hadronic medium without loosing much information. Thus, the lepton pair spec- 
troscopy in QCD will tell us a lot about the detailed structure of the hot/dense 
matter, which is quite similar to the soft-mode spectroscopy by the photon and 
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neutron scattering experiments in solid state physics. 
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Figure Captions 

Fig.l The light quark condensates in N=Z nuclear matter in the linear density 
approximation. Theoretical uncertainty of the nN sigma term is taken 
into account. We take y = 0.12 for the OZI breaking parameter, where 
y = 2{ss) n I (uu + dct)N with (-) N being the nucleon matrix element. 

Fig. 2 Borel curve for the p{ui) meson mass. Solid, dashed and dash-dotted lines 
correspond to p/p = 0, 1.0 and 2.0 respectively. Sq(p) determined by the 
Borel stability method at each density is also shown in GeV 2 unit. The 
Borel window is chosen to be 0.41GeV 2 < M 2 < 1.30GeV 2 . 

Fig. 3 Masses of p, u and 4> mesons in nuclear matter predicted in the QCD sum 
rules. The hatched region shows theoretical uncertainty. 

Fig. 4 Feynman diagram of the relativistic Hartree approximation (RHA). The 
wavy (dashed) line is the vector (scalar) meson propagator. The double 
(single) solid line denotes the Hartree (free) propagator of the nucleon. 

Fig. 5 An illustration of the divergent parts of E s . The dashed line is the scaler 
meson propagator and the double (single) solid line denotes the Hartree 
(free) propagator of the nucleon. 

Fig. 6 An illustration of the contribution from counter terms. Crosses denote 
the insertion of the counter terms with coefficients «i — 014. 

Fig. 7 (a) The Dyson equation for the vector meson propagator with self-energy 
n. The double (single) wavy line is the full (free) vector-meson propagator, 
(b) Vector-meson self-energy in the one-loop approximation. The double 
solid line denotes the nucleon propagator in the Hartree approximation. 

Fig. 8 Effective uj meson mass as a function of the baryon density using three 
different subtraction procedures for meson self-energy. The solid, dashed 
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and dash-dotted lines are obtained by the subtraction schemes given in 
ref.[P9"|,P0|j and [P^j , respectively. Triangle, diamond and square denote the 



values at p = p i n each subtraction scheme. 

Fig. 9 Effective masses of the nucleon and the uo meson as a function of the 
baryon density. m* e (m) denotes the real mass in the medium (the mass in 
vacuum). The dashed line corresponds to the invariant mass in the medium. 

Fig. 10 Real and invariant masses of the p-meson in the parameter set I, (g p , k p ) = (2.63, 
The dashed line corresponds to m* inv /m. 

Fig. 11 Same quantities with Fig. 10 in the parameter set II, (g p , k p ) = (2.72,3.7). 
The dashed line corresponds to m* nv /m. 

Fig. 12 Ratio of the baryon mass in matter M B and that in vacuum Mb (B = 



AT, A, S) as a function of p/po- This figure is taken from ref. ||31|[. 

Fig. 13 Ratio of the 0-meson mass in matter m*. and that in the vacuum as 
a function of p/po- Only hyperon contributions are included in the 0-meson 
self-energy. This figure is taken from ref. ||3T | . 



Fig. 14 Same with Fig. 13 except that only the nucleon contribution is included. 



This figure is taken from ref. |3T| 
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